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Abstract 

A long-standing issuei - ^ in the physics of strongly correlated electronic systems is whether the 
motion of a single hole in quantum antiferromagnets can be understood in terms of the quasiparticle 
picture. Very recently, investigations of this issue have been within the experimental reach^. 
Here we perform a large-scale density matrix renormalization group study, and provide the first 
unambiguous numerical evidence showing that in ladder systems, a single hole doped in the Mott 
antiferromagnet does not behave as a quasiparticle. Specifically, the injected hole is found to 
be always localized as long as the leg number is larger than one, with a vanishing quasiparticle 
weight and a localization length monotonically decreasing with the leg number. In addition, the 
single hole self-localization is insensitive to the parity (even-odd) of the leg number. Our findings 
may advance conceptual developments in different fields of condensed matter physics. First of 
all, the intriguing self-localization phenomenon is of pure strong correlation origin free of extrinsic 
disorders. Therefore, it is in sharp contrast to the well-known Anderson localization 17 and recently 
found many-body localization 18 , where extrinsic disordered potentials play crucial roles. Second, 
they confirm the analytical prediction a 14 ' 15 of the so-called phase string theor y 13 ' 14 , suggesting that 
the phase string effect lies in the core of the physics of doped Mott antiferromagnets. 
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The doped Mott insulator is generally believed to be the prototypical model of many 
realistic strongly correlated systems notably high T c cuprates^^i, and is a fundamental 
problem of great challenge in condensed matter physics. Important insights into doped 
Mott insulators may be gained by inspecting a special case namely the single hole doped 
Mott insulator. This has been under intense investigations^— for over decades, but re- 
mained highly controversial. On the one hand, a quasiparticle weight was found by using 
either analytic approaches (in combination with the self-consistent Born approximation)^~- 
or finite-size exact diagonalization^ on lattices up to 32 sites^. This indicates that an injected 
hole with a Fermi momentum maintains coherent motion in the quantum antiferromagnetic 
spin background. On the other hand, the validity of such a quasiparticle picture has been 
seriously questioned^ 1 ^— by various authors. It was first argued by Anderson 2 that the 
quasiparticle weight vanishes when a hole is injected into the Mott insulator. Later, this 
crucial observation received justification on a completely microscopic level: by identifying 
the so-called phase string effect^^, it was proved that the quasiparticle weight indeed van- 
ishes at the ground states. In addition, it was further predicted^ that the non-perturbative 
phase string effect may lead to self-localization of the injected hole. 

On the experimental side, in earlier angle-resolved photoemission spectroscopy (ARPES) 
studies^"—, broad spectral features have been observed in materials such as Ca2Cu02Cl2 
and Sr2Cu02Cl2- This shows indirectly that the motion of a single hole in quantum antifer- 
romagnets may not be understood in terms of the quasiparticle picture^ 1 ^. Very recently, 
an unprecedented degree of control has been reached in experiments with the Ca2Cu02Cl2 
parent Mott insulators, opening a new route for experimental study of a single charge doped 
Mott insulator. In particular, scanning tunneling microscope (STM) experiments on the 
atomic scale electronic structure of this material have unveiled a striking phenomenon: a 
single electron donated by surface defect creates an electronic state strongly localized in 
spaced. Consistent with the earlier ARPES experiments above, the observed localization of 
the single doped charge also does not support the quasiparticle picture. 

While numerical simulation serves as a powerful tool for resolving this puzzle, a crucial 
point is how to extend results obtained from finite-size numerical simulations to the large 
sample limit. To the best of our knowledge, this issue has not yet been addressed in the 
literature. The purpose of this work is to carry out the first large-scale numerical study and 
resolve the important issue of the validity of the quasiparticle picture in lightly doped Mott 
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ant if erro magnet s . 

To this end, we adopt the unbiased density matrix renormalization group (DMRG) 
method^ to extensively investigate this problem in ladder systems. We find that if a ladder 
sample is sufficiently long, the injected hole is always localized in the quantum antiferro- 
magnetic spin background, once the leg number is larger than one. In fact, the localization 
length gets continuously reduced with the increase of the leg number, pointing to a strong 
localization in the two-dimensional limit. Contrary to this, if the sample is sufficiently short, 
the injected hole remains itinerant. This suggests that quasiparticle behavior found in earlier 
numerical study 8,9 is likely a small-size effect. The intriguing phenomenon of self-localization 
of a single hole in Mott antiferromagnets and the vanishing quasiparticle spectral weight in 
the large sample limit constitute the unambiguous numerical evidence showing that a single 
hole doped into the Mott antiferromagnet no longer behaves like a quasiparticle. These 
support the analytic predictions^ 1 ^ of the phase string theory for doped quantum antifer- 
romagnets. 

Model and numerical method. — The DMRG calculation is performed for the t-J model 
Hamiltonian, 

H t j = ~t ( c i* c j* + h.c.) + J Y,( S i ■ S j ~ \ n i n i)- (!) 

(ij)" (ij) 

Here, c ia is the electron creation operator at site i, Si the spin operator, and n, the number 
operator. The summation is over all the nearest-neighbors, (ij). The Hilbert space is 
constrained by the no-double-occupancy condition, i.e., n» < 1. At half-filling, m = 1, the 
system reduces to Mott insulators (antiferromagnets) with a superexchange coupling, J. 
Upon doping a hole into this system, J2i n i — N — 1 (N the number of the lattice sites), and 
the hopping process is triggered as described by the first term of equation (JTJ, with t the 
hopping integral. 

Below, we shall focus on bipartite lattices of N = N x x N y , where N x and N y are the 
numbers of sites along the x and y directions, respectively. We shall study the ladders with 
a finite N y (from 1 to 5) and sufficiently large N x . We set J = 1 as the unit of energy and 
focus on the tj J — 3 case unless otherwise specifically stated. For the numerical simulation, 
we use fully open boundary conditions, and keep enough number of states in each DMRG 
block. Excellent convergence is achieved with total truncation error of the order of < 10~ 7 . 

Self-localization of the doped hole. — One of the main findings of the present work is that 
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the single hole is found to be localized for the leg number N y > 1. Examples of the hole 
density distribution, (n^) = 1 — (rij}, are shown in Fig. [1] (a) and (b) for N y — 3 and 4, 
respectively, where the sample size N x = 200 is clearly much larger than the localization 
lengths. In Fig. [TJ (a) and (b), we plot the hole density at a middle leg of the ladders along 
the x direction. We have checked that upon summing up the distribution from different 
legs, the sum rule J2i i^i) = 1 is satisfied. Examples of the contour plot of (n^ in the x-y 
plane can be found in Supplementary Information. 

To determine the localization length, we note that each density profile is well fitted by 
a Gaussian distribution function. Then, the localization length is defined as the full width 
at half-maximum (FWHM) of the distribution. If the doped hole is itinerant, the FWHM 
should increase monotonically with the sample length (N x ); in contrast, if the doped hole 
is localized, the FWHM should saturate at certain sample length, with the saturation value 
giving the localization length. Fig. [TJ (c) shows that the FWHM for each of the ladders, from 
N y = 1 to 5, increases linearly for small sample lengths, while saturates for sufficiently large 
sample lengths (except the N y = 1 shown in the inset of Fig. [TJ (c)). These results 

indicate that the single hole is well localized in space with a Gaussian profile of the density 
distribution. In particular, Fig. [1] (c) clearly shows that the saturated FWHM at N y > 1 
monotonically decreases with the increase of the leg number N y , implying the localization be 
even stronger in the two-dimensional limit. By contrast, there is no indication of saturation 
in FWHM for long one-dimensional (N y = 1) chains with N x up to 300 (see the inset of 
Fig. 11(c)), which is consistent with the fact that the doped hole in one dimensions is known 
to follow the Luttinger liquid behavior instead of being localized^ 1 ^. 

Figure [1] (a) and (b) show that the hole is localized in the central region of the sample. 
This is due to the fact that an open boundary condition is used for the DMRG calculation. 
We also add an artificial local chemical potential //(= 0.08) at the site x = 31, and find that 
the hole density profile is moved and peaks at the impurity accordingly (see Fig. [TJ (d)). 
This suggests the robustness of localization. 

Although the self-localization of the doped hole is insensitive to the parity (even-odd) of 
the leg number, as shown in Fig. [1] (a) and (b) as well as in Supplementary Information, the 
hole distribution function (n^ does clearly exhibit a parity effect: for the even-leg ladders 
(N y = 2, 4), there are always small site-dependent oscillations on top of the Gaussian density 
profile of (nf V while they are absent for the odd-leg ladders (N y = 3, 5). We point out 
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Fig. 1: Self- localization of the doped hole. In the 3-leg ladder system of size iV=200 x 3 (a) 
and 4- leg ladder of iV=200 x 4 (b), single hole doped into the antiferromagnets is well localized 
as shown by the density distribution (nf) in a middle leg (along the x axis), (nf) is well fitted 
by the Gaussian function with the full width at half-maximum (FWHM) defining the localization 
length. The behavior of the FWHM depends on the width N y (c): for the ladders of N y > 1, the 
FWHM first increases linearly at small N x and then saturates at large N x , with the saturation 
values much smaller than the maximal N x (= 200) and monotonically decreasing with the leg 
number; for the one-dimensional (N y = 1) chain, the FWHM increases monotonically with the 
sample length without saturation (the inset of (c)). With the open boundary condition, the hole 
is naturally localized at the central region of the sample. Nevertheless, it can be easily shuffled by 
adding a small local chemical potential [i as shown in (d) (where fi = 0.08 at x = 31). 

that this can be attributed to the underlying distinct spin-spin correlations for the odd- and 
even-leg ladders at half-filling that follow a power-law (reflecting the absence of spin gap) 
and exponential (reflecting the presence of spin gap) decay, respectively (see Supplementary 
Information). In fact, with the increase of the hopping integral t, the spin- gap effect will 
be gradually reduced such that the parity effect eventually disappears at large t/J limit, 
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FIG. 2: t/J-dependence of the hole density distribution. For different t/J ratios, we plot 
{n^) at the middle chain along the x axis for the ladders of size N = 80 x 3 (a) and N = 80 x 4 (b). 
The localization length (the FWHM) is reduced with the increase oft/J (see the insets). Moreover, 
the oscillation in the 4- leg ladder system becomes smaller with the increase of t/J, indicating the 
convergence of the even- and odd- leg ladders. 

where the aforementioned oscillation in (n^ for the even- leg ladders is also diminished as 
illustrated in Fig. [2] (b). The localization length is actually monotonically reduced as the 
ratio t/J increases as shown in the insets of Fig. |2J It thereby suggests that the detailed spin 
dynamic behavior, governed by the superexchange J, is not essential to the hole localization. 

Vanishing quasiparticle spectral weight. — One prominent issue regarding the validity of 
quasiparticle picture concerns whether the quasiparticle spectral weight vanishes at the 
Fermi momentum. Despite of great efforts^ - -^ 1 ^, this remains highly controversial in the 
literature. In this part, we turn to study this subject by invoking the state of the art 
numerical method. The spectral weight is defined as 

Z k = \ (lpl-hole\ C k a \lp0-hole)\ 2 , (2) 

where \ipi-hoie) is an eigenstate of H t _j of momentum — k in the one hole case, \ipo-hoie) 
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FIG. 3: Quasiparticle spectral weight. For both 3-leg (a) and 4-leg (b) ladders, the momentum 
distribution of the hole exhibits scaling behavior: after the rescaling, 1 — n(k) — > [1 — n(k)]N, 
the curves at different sample lengths in the insets collapse into a single one represented in the 
corresponding main panels. Note that for the 3-leg case, we fix k y = ^ and for the 4-leg case, 
k y = ^ and plot the distribution along k x . Representative contour plots in the whole k x -k y plane 
are shown in Supplementary Information. 

the ground state at half-filling, and Cfc CT is the Fourier transformation of the operator c, ia . 
Due to the translational symmetry of H t _j, is well-defined and characterizes the low- 
lying quasiparticle-like excitations. A finite Zy if means a finite overlap of the bare-hole state 
Cfc/o- \ipo-hoie) with the true ground state of a single hole at the Fermi momentum — kf. If 
Zk f — 0, then each injected hole will cause a global adjustment in the ground state, rendering 
the breakdown of the quasiparticle picture that is perturbatively tractable^. 

Importantly, the finding of the self-localization of the single hole implies that the doped 
hole should not behave like a conventional quasiparticle. Consistent with this observation, 
the quasiparticle spectral weight is numerically found to vanish, as we will show below. Tech- 
nically, different bases after truncation in the DMRG calculation of \ipo-hoie) and \ip\-hoie) 
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make it difficult to directly compute Zk f by using the definition (j2J). To overcome this 
difficulty, we note that a finite Zk f implies a sudden jump in the momentum distribution 
function n(fe) = J2a( c ka c kcr) at the Fermi momentum. Therefore, we first calculate n(k), 
that is a task well within the reach of the DMRG method, and then proceed to find Z^ r 
The results are shown in Fig. El 

The insets of Fig. [3] present the hole momentum distribution 1 — n(k) as a function of k x 
for fixed k y = 2ir/3 in the 3- leg ladder (a) and k y = n/2 in the 4-leg ladder (b). The value of 
k y is chosen in the way that the "sudden change" in 1 — n(k) can reach maxima by varying 
k x , according to the contour plots in the k x -k y plane (see Supplemental Information). It is 
very important that after the rescaling: 1 — n(k) — > [1 — n(k)]N, all the curves in the inset 
of Fig. E^a) [or (b)] collapse into a universal curve shown in the corresponding main panel. 
If one defines the Fermi surface by the sudden jump in the momentum distribution function, 
then the two universal curves suggest that the quasiparticle weight Zk f scales as 1/N for 
large N, and vanishes in the thermodynamic limit. 

To satisfy the sum rule, J2k[^ ~ n(k)] = 1, the width of the jump in 1 — n(k) must 
remain finite in the limit N — > oo. This is clearly shown in Fig. |3j consistent with a 
finite localization length in the real space. Such a finite width in the momentum space may 
also look like a remnant Fermi pocket of the hole (cf. the contour plots in Supplementary 
Information), whose density is finite in the localization volume even in the limit iV — > oo. 
We have also calculated the momentum distribution of the hole at different ratios of t/J. 
As shown in Supplementary Information, for a given sample size, the jump near the Fermi 
point gets continuously reduced with increasing t/J, which is consistent with earlier work 8 . 

Moreover, it has been predicted analytically^ that a spin-charge separation occurs as 
a spinless ho Ion is localized while a neutral spinon of 5* = 1/2 moves away in the two- 
dimensional case. The present numerical simulations indeed confirm this prediction for the 
odd-leg ladders where the spin gap vanishes. Specifically, we inject a hole into a 3-leg ladder 
by removing a j. spin electron. The results are shown in FigJU the hole (charge) is localized 
at the sample center whereas the extra spin of S z = 1/2 is spread along the x direction, 
with a (coarse grained) density (Sf) c approximately uniform. (For the even-leg ladders, 
the spin gap renders the observation of the spin-charge separation more difficult. Indeed, it 
is required that the spin gap to be small enough, a condition that can be achieved only if 
the leg number is sufficiently large.) 
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FIG. 4: Spin-charge separation. We remove a J, spin electron in a N = 200 x 3 ladder and 
observe spin-charge separation. The hole (charge) is localized at the sample center, with the density 
distribution of Gaussian type; the spin of S z = 1/2 spreads over the entire sample, with the (coarse 
grained) density (Sf} c . g . approximately uniform. Note that to compute the spin density, a local 
coarse graining has been performed so as to average out local antiferromagnetic oscillations. 

Discussion. — We note that in Fig. [TJ the localization lengths are much larger than 
the lattice constant. This indicates that the injected hole moves "freely" from one site to 
the other within the localization volume. At larger (than the localization length) scales, 
such a free motion is fully suppressed, leading to the hole localization. In this sense, the 
scenario resembles the quasi-one-dimensional Anderson localization discovered by Efetov and 
Larkin long time ago™ There, quantum diffusion occurring at short scales is brought to a 
halt by wave interference at large scales. The phenomenological analogy notwithstanding, 
the two systems are fundamentally different: in the present system, there are no extrinsic 
disorders at all. Instead, as we will show below, the self-localization of the injected hole 
in the quantum antiferromagnetic background arises from intrinsic disorders of pure strong 
correlation origin. Therefore, it is distinctly different from either conventional Anderson 
localization 1 - or many-body localization^ discovered more recently. 

Let us consider the motion of the hole in the quantum antiferromagnetic spin background 
from the injection point to a distant site (Fig. [5]). The wave nature allows the hole to 
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14 given a hole path, 



propagate along all the possible paths. As discovered in Refs. 
p, the quantum phase is completely determined by the parity of the hole- J, spin exchange 
number, N^lp}. More precisely, upon hopping to its nearest neighbor, the hole acquires 
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FIG. 5: Interference picture of self-localization. A sign ± is acquired when the hole exchanges 
its position with the nearest neighbor, an f (\.) spin. Consequently, as the hole moves from the 
injection point to a distant site, a sign sequence - the phase string - is left behind as exemplified 
by the green (solid) and dashed (blue) lines. These phase strings destructively interfere with each 
other, suppressing the forward scattering and leading to the self-localization of the injected hole in 
the quantum antiferromagnetic spin background. 

a sign, ±, depending on whether it exchanges its position with f or ^ spin. As the hole 
moves along a given path, a sign sequence known as the phase stringi^ 4 is left behind. The 
quantum amplitude characterizing the hole motion is the superposition of all the propagation 
amplitudes associated with individual hole paths, and is given by 

EEfr^}](-i) <bl . (3) 

P {</,} 

where the first sum is over all the possible hole paths, p, the second over all the intermediate 
spin configurations, {0}, encountered in the path, p, and p P [{0}] > 0. Equation ([3]) suggests 
that as a unique property of Mott physics, the spin configurations provide a translationally 
variant background - the analog of disorder potentials - for the hole motion. 

To better understand the "disordering" nature of the phase strings, we note that as the 
hole path is slightly distorted, the intermediate spin configurations encountered may also be 
slightly changed. But the resulting change in the parity of the number can be signifi- 
cant. That is, the "non-integrable quantum phase factor", (—1)^^, will strongly fluctuate 
as the hole paths are sampled. Then, from equation ([3]), we see that even paths close to 
each other can destructively interfere very singularly. As a result, the forward scattering 
is significantly suppressed, and eventually, localization is reached when the system is suffi- 
ciently large. Indeed, as the leg number or the parameter t/J increases, such destructive 
interference proliferates, resulting in the decrease of the localization length (Fig. [1] and [2]) . 
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In contrast, provided that the leg number is unity, the destructive interference from differ- 
ent paths no longer exists and the hole escapes from localization, in agreement with the 
numerical simulation (inset of Fig. [I]). For the leg number comparable to the sample length, 
the destructive interference arising from the disordered signs leads to two-dimensional self- 
localization, as proved analytically in Ref. |l5|. Therefore, we conclude that the intrinsic 
"disordered" sign structure, (—1)^^, is responsible for the self-localization of the injected 
hole and the vanishing quasiparticle spectral weight. 

In summary, we provide the first unambiguous numerical evidence showing that a sin- 
gle hole injected into Mott antiferromagnets is localized in the quantum antiferromagnetic 
background, and that the quasiparticle weight vanishes in the thermodynamic limit. These 
findings are consistent with both the analytic prediction of the phase string theory for a sin- 
gle hole-doped two-dimensional Mott antiferromagnet and the experimental results achieved 
by using ARPES and STM methods. In view of that the phase string effect is a rigorous 
property of the t-J model on bipartite lattices regardless of doping concentration, temper- 
ature, and dimensions 2 ^, the present work may provide significant new insights into the 
long-standing issue of doped Mott insulators and thereby high temperature superconductiv- 
ity. 
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DMR 1121053), by the NBRPC (nos. 2009CB929402, 2010CB923003, 2011CBA00302 and 
2011CBA00108), and by the Tsinghua University Initiative Scientific Research Program (no. 
2011Z02151). 
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At half-filling, the t-J model reduces to the Heisenberg model. For the isotropic Heisen- 
berg coupled-chain systems, the behavior of the even-leg ladders is dramatically different 
from that of the odd ones. A well-known fact is that the even-leg ladders have a spin gap, 
while the odd-leg ladders are gapless, leading to an exponential decay of the spin-spin cor- 
relation in the former, while a power law decay in the latter (see Fig. 1ST]) . The spin gap for 
the even-leg ladders is expected to vanish in the large leg number (namely two-dimensional) 
limit. The spin structure factor for the 4- leg case has already shown strong antiferro mag- 
netic correlations as illustrated in Fig. [S2j These results are consistent with earlier DMRG 
work 1 . 
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Fig. SI: The main panels show the spin-spin correlation 
located on the middle leg of 3-leg (a) and 4-leg ladder (b), respectively. The insets are In- In plot in 
(a) and ln-linear plot (b), fitted with the straight (blue) lines which indicate the power law decay 
for the former and exponential decay for the latter. 
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Fig. S2: Spin structure factor for the 4-leg ladder (iV=80x4), in which the peak is located at (vr,7r) 
indicating strong spin antiferromagnetic correlations even though a finite spin gap is present. 

For the hole density distribution (n^) calculated in this work, the total hole number is 
fixed, J2{ n i) — 1- The contour plot (nf) in the x-y plane of the ladder is illustrated in 

i 

Fig. [S3] for N = 40 x N y with N y = 2, 3, 4, and 5. Besides the localization in the central 
region along the chain (x) direction, a prominent feature in the profiles is the distinction 
between the even- and odd-leg ladders: the former has a spatial oscillation, while the latter 
has none. This is concomitant with the parity effect in the spin-spin correlation seen in 
Fig.EU 

Fig. [S4] shows the contour plot of the electron momentum distribution function n(k) for 
the 3-leg (a) and 4-leg (b) cases, respectively at N x = 80. The minimum of n(k) appears 
at k y =±7i/2 for the 4-leg ladder and k y =±2n /3 for the 3-leg ladder. The jump of the hole 
momentum distribution as a function of the ratio t/J is shown in Fig. IS5l which indicates 
that the jump is monotonically reduced with the increase of t/J . 



1 S.R.White, R. M. Noack, and D. J. Scalapino. Phys. Rev. Lett. 73, 886(1994). 
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Fig. S3: Contour plots of the hole distribution in real space. The sample size iV is 40 x 2 (a), 40 x 3 
(b), 40 x 4 (c), and 40 x 5 (d) from top down. 




(a)N=80x3 (b)N=80x4 

Fig. S4: Contour plots of the electron momentum distribution n(k) for 3-leg (a) and 4-leg (b) 
ladders. 
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Fig. S5: The hole momentum distribution at different ratios of t/J for 3-leg (a) and 4- leg (b) 
ladders. 
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